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(57) AbrSge/Abstract 

A method and device for valuation of financial derivatives, wherein a value of a derivative is computed by a determination of an 
expectation. Input parameters are communicated by an input unit to a computer, such as at least one processor, to establish an 
integrand as a function of the input parameters. A multivariate integration domain is computed. A sparse grid method is used to 
determine integration points and integration weights as a function of the input parameters. The integrand is integrated with an 
integration domain at the integration points to determine integrand values. One or more expectation parameters are computed 
by combining the integrand values and the integration weights. A value of the derivative is communicated through an output unit, 
for example to a display monitor or another display device. 
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ABSTRACT OF THE DISCLOSURE 
A method and device for valuation of financial derivatives, wherein a 
value of a derivative is computed by a determination of an expectation. Input 
parameters are communicated by an input unit to a computer, such as at least one 
processor, to establish an integrand as a function of the input parameters. A 
multivariate integration domain is computed A sparse grid method is used to 
determine integration points and integration weights as a function of the input 
parameters. The integrand is integrated with an integration domain at the integration 
points to determine integrand values. One or more expectation parameters are 
computed by combining the integrand values and the integration weights. A value of 
the derivative is communicated through an output unit, for example to a display 
monitor or another display device. 
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METHOD A ND DEVICE FOR EVALUATION OF FINANCIAL 
DERIVATIVES USING SPARSE GRIDS 

BACKGROUND OF THE INVENTION 
Field of the Invention 

This invention relates to a method and device for electronic data 
processing, particularly in the field of financial informatics. More particularly, this 
invention relates to fast and accurate valuation of financial derivatives using 
electronic computer systems. 

Description of Related Art 
The valuation of financial derivatives has become popular in die last 
decade and has become an important instrument in financial economics. The 
possibility of such valuation is a necessary prerequisite for a financial institution to 
be able to offer a new derivative. Also, this valuation is necessary in order to 
structure, collateralize, and optimize a financial portfolio. 

There are many different kinds of financial derivatives, such as those 
based on interest rates and/or assets. This includes, among others, bonds, swaps, 
future, CMOs, and options. Their valuation occurs under die assumption of the 
arbitrage principle using partial differential equations or the martingale approach. 
The martingale is more general, where a stochastic process for an underlying value, 
such as an interest rate or asset price, is first specified. Then, the equivalent 
martingale measure which converts the underlying process into a martingale is 
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determined Finally, the value of the derivative is computed as the expectation of its 
discounted payoff function under this risk-neutral measure* Under certain conditions 
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the approaches via partial differential equations and martingales are equivalent The 
connection is made by generalized Feynman-Kacz formulas. However, the 
martingale approach is more universal and easier to adapt to new situations. 

The price of a financial derivative can thus be expressed as an 
expectation. In continuous time, the integrand itself contains a path integral, which 
can be discretized with the Euler method/trapezoidal rule or similar methods. 
Sometimes, induced by the application, a time-discrete model is used Both cases 
lead to high-dimensional integration problems. For the simple European call option 
it is possible to give a closed-form solution, a Black-Scholes formula, but more 
complex options require a numerical solution method This holds for other types of 
financial derivatives analogously. 

Classical multivariate quadrature is not suited as a numerical integration 
method for high-dimensional integrands. Problems associated with dimension are 
encountered because the work scales exponentially with the dimension. The 
complexity is of the order O (N***) where r is the smoothness of the integrand and d 
its dimension. On the other hand, the Monte Carlo method is independent of the 
dimension. Here, die integrand /is evaluated at a random series of N points Xj which 
results in the following quadrature formula. 

Q»f = )j\Ax\ Equation 1 
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The Monte Carlo methods converges very slowly and only in a 
stochastic sense. The accuracy which can be reached with # function evaluations is 
of the order 0(1/V5y). With Quasi-Monte Carlo methods which have been developed 
in the last decade, the integrand is evaluated at a deterministically determined series 
of points Xj and in analogy to the Monte Carlo method, the following quadrature 
formula is used 

QX-^i/h) Equation2 

There is a variety of different constructions, for example the Halton, 
10 Sobol sequences or Faure sequences, which differ in corresponding pre-asymptotic 

behavior and all have an order of convergence of O ((log Nf/N). In addition the error 
is deterministic. Prototype is the program FinDer of J. Traub, S. Paskov, Faster 
evaluation of financial derivatives. Journal of Portfolio Management 22, 1, 113-120 
(1995), which is used by many banks. This method is disclosed by United States 
Patents 5,940,810 and 6,058,377. 

On the other hand, the so-called sparse grid method is an approach 
where multivariate quadrature formulas are constructed by suitable combination of 
tensor products of univariate quadrature formulas, such as the Clenshaw-Curtis or 
Gauss-Patterson formulas. 
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The general sparse grid metb^ Consider 
a series of univariate quadrature formulas for a univariate function/ represented by 
the following formula. 



0!/:«2>.-/tr.) Equation3 

Now, define the difference formula by 

Ai/ := {Ql-QL)f with Eqa^onA 
Qj/ := 0 

The sparse grid construction for ^-dimensional functions f consists for 
7eINandkeIN < 'of 



Q?f :=y L{*\ • .«» A*)/ Equations 



with index sets // such that die following function holds for all kel f 

k - e y € /, for Isys d, kj> 1 Equation 6 

Special cases of the method are classical sparse grids, where l t = {|k| , 
zl + d- 1 }, as well as classical product formulas where 7 /== { | k . Fig.4shows 
examples of three different classical sparse grids in the 2D case based on the 
trapezoidal rule, the Clenshaw-Curtis formula, die Gauss-Patterson formula, and the 
Gauss-Legendre formula. 

One substantial difference to Monte Carlo and Quasi-Monte Carlo 
methods is the use of weight factors of different size. Now, in die representation 
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Q*f~ 2* /«= i w if\ x * ) p the weights w t of the sparse grid method are not all equal to 
UN. 

The order of convergence of the classical sparse grid method is 
8 = O0°g(N) (d " lXr4 * l) N" 1 ) and thus it is also independent of the dimension d. In contrast 
to Monte Carlo and Quasi-Monte Carlo methods, the sparse grid method is able to use 
the smoothness r of the integrand and thus has exponential convergence for smooth 
integrands (r - «>). Thus, for smooth integrands this method is substantially faster 
than the Monte Carlo or Quasi-Monte Carlo methods. 

For the sparse grid method, one problem is that die convergence rate 
deteriorates for non-smooth integrands which show up for options, and its advantage 
is thus lost Another problem is that the method is, just as the Quasi-Monte Carlo 
method, not entirely independent of the dimension, and die rate of convergence also 
degrades with rising dimension. 

SUMMARY; OF THE INVENTION 

It is one object of this invention to provide a method and a device for 
fast and accurate valuation of complex financial derivatives, such as for use by banks, 
insurance companies and/or private investors. 

It is another object of this invention to provide a method and device for 
determining a price of the complex financial derivatives, as a function of an 
expectation. 
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It is yet another object of this invention to provide a method and device 
that increases accuracy and decreases running time as required by conventional 
devices and methods for computing the expectations, such as based on Monte Carlo 
or Quasi-Monte Carlo integration. 

The above and other objects of this invention arc accomplished, 
according to one embodiment of this invention, with a method and device, such as a 
computer, that performs valuation of financial derivatives based on multivariate 
numerical integration using sparse grids. The method of this invention can compute 
arising expectations more efficiently than conventional methods and thus can compute 
and determine prices of the derivatives much faster and with much more accuracy 
than conventional methods and devices. 

In one embodiment of this invention, the device includes an input unit, 
an output unit, a memory unit and a program unit The program unit has a setup 
module, a discretization module and an integration module. In die event of non- 
smooth integrands, for example those associated with option pricing, transformations 
decompose die integration domain into smooth subdomains. Dimension reduction 
and parallelization are also used to accelerate the method of this invention. 

The method and device of this invention provides fast and accurate 
determination of objective prices and subjective prices of financial derevatives, such 
as bonds, swaps, futures, CMOs, and options. The calculated prices can be used in 
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purchasing or selling the derivatives, such as for portfolio optimization, pricing of 
contracts, risk management and/or investment decision making. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Features of the method and device according to this invention can be 
better understood when die specification is read in view of the drawings, wherein? 

Fig. 1 is a block diagram showing different components of a device 
according to one embodiment of this invention; 

Fig. 2 is a block diagram showing a subdomain, processed in a parallel 
computer system with distributed memory, according to one embodiment of this 
invention; 

Fig. 3 is a block diagram of a subdomain processed in a parallel 
computer system with shared memory, according to another embodiment of this 
invention; 

Fig. 4 is a diagram showing two-dimensional renderings of classical 
sparse grids based on the trapezoidal rule, the Genshaw-Curtis formula, the Gauss- 
Patterson formula, and the Gauss-Legender formula; 

Fig. 5 shows diagrammatic views of three-dimensional renderings 
which compare different integration methods for valuation of path-dependent options; 
and 

Fig. 6 shows graphic representations of an error parameter versus a time 
parameter for different methods, 

8 
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DESCRIPTION OF PREFERRED EMBODIMENTS 

This invention relates to a device comprising a computer and a method 
which is based on the sparse grid method This invention allows the fast computation 
of expectations which are required for the valuation of derivatives from finance and 
banking, for example, interest rate, asset, currency, and/or commodity derivatives. 

Fig. 1« shows a block diagram of die method and device which 
comprises a computer having an input unit (1) and an output unit (6), a memory unit 
(5) and a program unit (2-4). The input and output units (1, 2) facilitate the input of 
the parameters of the derivative and the output of the value of the derivative. The 
memory unit is used for storing the program, die input and output values, and 
intermediate results. The program unit (2-4) comprises a setup module (2), a 
discretization module (3), and an integration module (4). The setup module (2) is 
used to determine the integrand function. The discretization module (3) determines 
suited integration points using sparse grids. The integration module (4) computes the 
value of the derivative through combination of the function values at the integration 
points. 

The following specification more precisely describes the individual 
components. The input (1) can be fed by any suitable input, such as a keyboard, 
another input device, or an analog or digital connection. The input data specifies the 
exact type of derivative as well as the properties of underlying values. This data as 
well as other intermediate results are stored in memory. A computer program reads 
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this data in a first step, refered to as setup and constructs a multivariate integration 
problem based on the input data. In a second step, referred to as discretization, suited 
integration points and integration weights are determined based on the sparse grid 
method using this information. In die third step, refered to as integration, the 
integration problem is solved through evaluation of an integrand at computed 
integration points and by combination of the results. The value of the derivative is 
returned as the final result The output (6) can be read through a monitor, another 
output device, or an analog or digital connection. 

The method of this invention can be used to smooth integrands which 
10 typically arise in the valuation of interest rate derivatives without further 
modification. Examples are portfolio optimization, the pricing of futures and forward 
contracts, the valuation of mortgage-based securities, the pricing of insurance 
contracts, risk management, or investment decision making. In the case of non- 
smooth integrands, for example which arise in option pricing, the convergence rate 
of the method is deteriorated. In practically all cases the integrands have 
discontinuous derivatives (r — 1), in some cases the integrand itself is discontinuous 
(r = 0). According to this invention, there is a decomposition of the integration 
domain into smooth parts in the discretization step;- for example by zero finding. 
Using suited transformations the sparse grid quadrature formulas are mapped to these 
20 subdomains and the total integral is computed as the sum of these partial integrals. 
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Thus, integration takes place only over smooth areas and the advantageous properties 
of the sparse grid method are maintained 

The sparse grid method is.largely, but not totally, independent of the 
dimension of the integration problem. The dimension d occurs in the order of 
convergence as the exponent of a logarithmic factor. This leads to a relatively slow 
degradation of the convergence as the dimension d is rising. Therefore, it is 
reasonable to use dimension-reducing methods. In the case of path-dependent options 
the underlying stochastic process can be discretized hierarchically, for example using 
the Brownian bridge construction. In the case of performance-dependent options, die 
volatility matrix is transformed using singular value decomposition. In both cases, 
a reduction of die effective dimension is achieved by focusing on the respective most 
important dimensions. In die discretization module, die sparse grid method uses more 
accurate quadrature formulas in more important dimensions and less accurate 
formulas in less important dimensions. This way, in many cases a practical 
independence of the dimension is achieved for die method. In other cases at least a 
strong acceleration of the method can be obtained. " 

An additional acceleration of the method of this invention occurs 
^through parallelization, such as a distributed computation on several processors. At 
this step, the entire integration problem is decomposed into independent subproblems. 
These are, for example, the subdomains originating from above transformation or the 
partial sums over the different multiindices k in the sparse grid quadrature formula. 
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The device is shown in Figs. 2 and 3. The distribution module (7) performs the 
decomposition of the integration problem. The combination module (8) combines the 
partial results during the parallel computation. The individual partial problems are 
processed either in a parallel computer system with distributed memory, such as 
shown in Fig. 2 or with shared memory, such as shown in Fig. 3. 

EXAMPLE 

An option is a contract between two parties which gives the holder the 
right but not die obligation to buy/sell a fixed amount of assets at a certain date for a 
fixed price. Thus, an option has an intrinsic value, since the holder does not have to 
exercise the option. One question to answer is how to obtain a fair price for the 
option. More precisely, using the martingale approach there is a value: 

V=e- r E(p({sf})) Equation? 

Here, r ii the constant interest rate, / = 7* is the exercise time, P is the 
S J 

payoff function, and 1 is the asset price of the /-th asset of the portfolio at time jut 

with a/ = T/M. Here, M is the number of discrete time stype j — 1, and Af and 
each denotes the number of considered assets /= 1, „ y N m e^is the discount factor and 
/ = 0 is the time of valuation. 
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Assuming that the assets of die portfolio follow a model described by 
the stochastic differential equation 

N 

^ = M^ + X^A^ Equations 



where (a^) denotes the time-constant volatility matrix, the time-constant drift for 
asset / and dWt a geometric Brownian process for asset z, such as W { is the Wiener 
measure. There are also approaches with time-varying volatility matrices and time- 
varying drift In this case a* and # follow own SDEs modeling their behavior. 
Integration and Itd-formula yield after time discretization using the time steps 7 a/ the 
values sL 

Sf =Sf e Equation9 
and thus . - v . ^ 

^ y ± ^0 Y. Jl* -? II ^JJ */+X m a a JAJf^J Equation 10 

Here * are N{0 9 l)-distributed random variables. In a computer program the first 

recursive formula is used However, following non-recursive summated form is 
necessary to define the later integrand. 
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The definition of the expectation £(.) is now for a general fiinction / 

defined by: 

MMh £ - £ & \) .-g(z?)dzi ~,d X » Equation 11 



with the standard normal distribution 



= — L. e -?n Equation 12 



With the requirement that arbitrage does not exist, an expectation is 
taken with respect to the equivalent martingale measure of the underlying process. 
Here, fi ( is replaced by r in the integral of the expectation. This way, also the change 
from the stochastic differential equation into its risk-neutral form takes place. 

Now, only the concrete payoff functions are missing and these depend 
on the respective option type. In general, the functions have the structure 

P({S/}) = ma^O,^^/}) - k) Equation 13 

for call options and 

P({S/}) = ma^O,^-^/})) Equation 14 

for put options. Thereby, K is the strike price. Examples are path-dependent options 
where \ 9 M>\ and performance-dependent options where N>1, M= 1. 
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The most popular example for path-dependent options are so-called 
Asian options. Here, die average over the asset prices at all points in time between 
the current time and the exercise time are taken for the payoff function. Consider the 

case of the geometric average such as U My V # ' * ' i ) m Jaffas 

case, a closed-form solution in form of a generalized Black-Scholes formula is 
available. Small changes and variations, for example already the arithmetic average, 
are not solvable analytically and require a special numerical integration method 

An example for performance-dependent options are conventional so- 
called Tailored options". Here, the price of the option depends on the relative 
performance of one asset in comparison to the other assets in a considered set of 
assets. One example for a call option is 

a max(s' l - K,6) else 



where a is a parameter in [0, 1]. A mixture of both option types is, of course, also 
possible. M- N= 1 yields as a special case the European call option with the Black- 
Scholes formula as analytic solution. 

Now, for option pricing it is generally necessary to solve a high- 
dimensional integration problem. The dimension results from the time steps over 
which the average is taken for patb-dependent options and/or from the number of 

15 
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considered assets for perfonnance-dependent options. The integration domain is so 
far due to the definition of the expectation. In order to be able to apply 

integration methods, the integral of the expectation is transformed using the 

f 

cumulative normal distribution G(y) ~ g(x)dx onto [0,1^ to get 

444)= C- £ /(M4)*: - *i Equad00ls 

For the integration of Equation 16 with Equation 13 or Equation 15 the 
sparse grid method is applied profitably. 

For options one problem is that the payoff function is no longer smooth 
due to the nature of an option. This ultimately expresses that the option will not be 
exercised if the purchase or sale of the underlying asset would result in a loss. The 
integrand shows a kink, path-dependent options, or even jumps, performance- 
dependent options, with respect to a (Af • N r l)-dimensional manifold Besides the 
singularities on the boundary of the integration domain which result from the 
transformation to [0,1]* Examples for integrands in the two-dimensional case are 
shown in Fig. 5. The kink, respectively jump is clearly visible. The smoothness 
requirements for die application of the sparse grid method are no longer fulfilled. 

Now integrating only over the support of the integrand, the integrand 
is a smooth function, and the kinks and jumps are exactly located on the boundary of 
the support To determine the integration area it is sufficient to compute the zeros of 

16 
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the integrand When the integral is computed iterativeiy, the zero finding is restricted 

to one, the last, dimension. The zero % in the last dimension, Newton method for the 

kink or bisection for the jump, is determined and the integrand with respect to the last 
dimension is transformed with the linear mapping 

Az)=Z •(!-/)+ X Equation 17 

onto [0,1]. 

Fig. 5 shows a comparison of the different integration methods for the 
valuation of a path-dependent option with 6 times steps, M = 6, N « 1, and a 
performance-dependent option with M = 1, N= 2. The superior rate convergence of 
sparse grid method using transformation and Gauss-Patterson formulas with respect 
to other methods, Monte Carlo (MC), classical product approach without (PR) and 
with transformation onto the support (PRTR), Quasi-Monte Carlo without (QM) and 
with transformation onto the support (QMTR), as well as Gauss-Patterson sparse grids 
without (SG) and with transformation onto the support (SGTR), is clearly visible. 
The error is plotted logarithmically and shows that the method of this invention is 
superior to the Monte Carlo and Quasi-Monte Carlo methods which cannot use the 
higher smoothness of the integrand while the general sparse grid quadrature of this 
invention can use it in an optimal way. 

Since the complexity of the sparse grid method of this invention is not 
totally independent of the dimension of the problem due to the term log (N) (d ~ n{r * n , it 
makes sense to also apply dimension reduction methods. The method can be 

17 
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accelerated using a hierarchical discretization of die stochastic process and adaptive 
refinement 

The most natural way to discretize a stochastic process is through a 
random walk, such as by the recursive formula 

^y+t _ gj e *{ w **) Equation 18 

where exactly corresponds to the exponent in Equation 19. In the Brownian 

bridge discretization, the process is discretized using a fixture and past value 

, _ Sf* k + S{ * MSuw/) Equation 19 

* 2 



This way, starting with °* 9 °* z 



the values 



rrA//2 C A//4 rO* /4 r*3M/S 



etc. are determined. This leads to a 



concentration of die total variance of the process in the first steps of the discretization 
which improves the convergence rate of Quasi-Monte Carlo methods. 

For the classical sparse grid method there is no immediate advantage 
from this discretization technique since all dimensions are of same importance. The 
general sparse grid method can be applied dimension-adaptively, and quadrature 
formulas with lower degree can be used in less important dimensions. Thus, die 
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effective dimension of the integrate such as ^ 

derivatives and the resulting integration problems can be computed fester and more 
accurately. 
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WE CLAIM: 

1. In a method for valuation of financial derivatives, wherein a 
value of a derivative is computed by a determination of an expectation, die method 
comprising: 

inputting aphirality of input parameters of the derivative to at least one 
processor and establishing an integrand as a function of die input parameters; 

computing a multivariate integration domain; 

applying a sparse grid method to determine a plurality of integration 
points and a plurality of integration weights as a function of the input parameters; 

evaluating the integrand inside an integration domain at the integration 
points to determine a plurality of integrand values; 

computing an expectation by combining die integrand values and die 
integration weights and determining a value of the derivative from die expectation; 
and 

outputting the value of the derivative. 

- 2. In the method according to Claim 1 wherein the input parameters 
and die value are communicated as at least one of a digital signal and an analog 
signal. 
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3. In the method according to Claim 1 wherein the integration 
domain is determined by finding an integration area by computing zeros of the 
integrand. 

4. In the method according to Claim 1 wherein the integration 
points and die integration weights are determined dynamically or adaptively when 
evaluating the integrand. 

5. In the method according to Claim 1 wherein the at least one 
processor computes the integration domain as a plurality of partial integrals* 

10 

6. In the method according to Claim 1 wherein the at least one 
processor computes in parallel at least one of the integration points, die integration 
weights, die evaluation of the integrand and die combination of the integrand values. 

7* In die method according to Claim 1 wherein a dimension of 
integrals is reduced to accelerate the method. 
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8. In themethod according to Oaim 1 wherein at least one second 
expectation is computed, and die expectation and the at least one second expectation 
are combined 

9. In the method according to Claim 1 wherein one of a keyboard 
input and an input device is used to input the input parameters. 

10. In the method according to Claim 1 wherein the integrand is 
evaluated by integrating the integrand 

10 

11. In the method according to Claim 1 wherein the value of the 
derivative is communicated as at least one of a digital signal and an analog signal and 
the value is displayed on at least one of a monitor and an output device. 



22 



CA 02364920 2001-12-12 



12. A device for valuation of financial derivatives, wherein a value 
of a derivative is computed by a determination of an expectation, die device 
comprising? 

an input unit communicating a plurality of input parameters of the 
derivative to a computer; 

the computer co mprisin g- a setup module establishing an integrand as 
a function of the input parameters and computing a multivariate integration domain, 
a discretization module applying a sparse grid method to determine a plurality of 
integration points and a plurality of integration weights as a function of die input 
parameters, and an integration module evaluating the integrand inside an integration 
domain at the integration points to determine a plurality of integrand values and 
computing an expectation by combining the integrand values and the integration 
weights; and 

an output unit communicating a value of die derivative. 

13. In the device according to Claim 12- wherein the input unit 
c ommuni cates the input parameters and the output unit communicates fee value as at 
least one of a digital signal and an analog signal. u*:*« 
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14. In the device according to Claim 12 wherein die computer 
comprises at least one processor computing the integration domain as a plurality of 
partial integrals. 

15. In the device according to Claim 12 wherein the computer 
comprises ^plurality of processors computing in parallel at least one of die integration 
points, the integration weights, the evaluation of die integrand and die combination 
of the integrand values. 

10 16. In die device according to Claim 12 wherein die input unit 

comprises one of a keyboard input and an input device for inputting the input 
parameters. 

17. In the device according to Claim 12 wherein die computer 
communicates the value of the derivative as at least one of a digital signal and an 
analog signal. 

18. In the device according to Claim 12 wherein the output unit 
comprises at least one of a monitor and an output device. 
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